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We demonstrate that a giant spin Hall effect (SHE) can be induced by introducing a small amount 
of Bi impurities in Cu. Our analysis based on a new 3-dimensional finite element treatment of spin 
transport shows that the sign of the SHE induced by the Bi impurities is negative and its spin Hall 
(SH) angle amounts to —0.24. Such a negative large SH angle in CuBi alloys can be explained by 
applying the resonant scattering model proposed by Fert and Levy [Phys. Rev. Lett. 106, 157208 
(2011)] to 6p impurities. 

PACS numbers: 72.25.Ba, 72.25. Mk, 75.70.Cn, 75.75.-c 



Spintronic devices manipulating pure spin currents, 
flows of spin angular momentum without charge current, 
should play an important role in low energy consumption 
electronics of next generation. This explains the current 
interest for the spin Hall effect (SHE) which provides a 
purely electrical way to create spin currents without fer- 
romagnets and magnetic fields. The SHE, originally pre- 
dicted by Dyakonov and Perel [l[ , can be described as an 
accumulation of spins generated by an electric current on 
the edges of a nonmagnetic conductor Q. Its interest in 
spintronics comes from its application to convert charge 
into spin currents (or vice- versa). The spin Hall (SH) an- 
gle, characteristic of the conversion yield between charge 
and spin, reaches a few % in heavy metals with strong 
spin-orbit (SO) interactions such as Pt 

In addition to the intrinsic SHE of pure metals [H, Q , 
the skew scattering [l(J and the scattering with side- 
jump jllj on impurities with strong SO interactions can 
also give rise to the SHE. Such extrinsic mechanisms have 
already been studied in the SHE induced by nonmagnetic 
impurities in Cu [12| and in the anomalous Hall effect of 
ferromagnetic alloys 13( . One of the big advantages for 
the extrinsic SHE is that it allows controlling the SH an- 
gle by changing the combination of host and impurity 
metals as well as by tuning the impurity concentration. 
According to recent theoretical predictions 14 , HI , some 
combinations of noble metals and impurities could induce 
very large SH angles, for example in Cu or Ag doped with 
Bi. In this work we obtain a large SHE signal in Cu doped 
with a small amount (< 0.5%) of Bi. From our analysis 
based on a 3-dimensional (3D) finite element treatment 
of the spin transport equations , the SH angle can be 
estimated to be —0.24 at T = 10 K, which is an order 
of magnitude larger than that in pure metals such as Pt, 
Pd and almost twice larger than that announced recently 
for the (3 phase of Ta [H]. 
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FIG. 1: (Color online) (a) Schematic of the ISHE measure- 
ment. The ISHE in CuBi deflects spin- up and down electrons 
|e| (e is the charge of the electron) denoted by spheres with 
arrows to the same side. Other arrows indicate the electron 
motion direction, (b) Electrochemical potential distributions 
for spin-up and down electrons in Cu. Spin-up and down 
electrons diffuse in opposite directions, which characterizes a 
pure spin current (no charge current), (c) Scanning electron 
microscopy image of our spin Hall device. The yellow arrow 
in (c) represents the positive direction of the magnetic field 
in the case of ISHE measurement. 
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Figure 1(a) shows the 
spin absorption method | 
trie current flows from Pyl to the left side of the Cu 
wire, the resulting spin accumulation induces a pure spin 
current (no net charge current, i.e. 7c = or If = — 1± 
for opposite flows of spin-up and down electrons) on the 
right side of the Cu wire [see Fig. 1(b)]. As discussed later 
on, most of the pure spin current is absorbed vertically 
into the middle wire below Cu. The opposite spin-up 
and down electrons composing the absorbed pure spin 
current are deflected to the same direction (along the x- 
axis) by the ISHE in CuBi, and an ISHE voltage Vishe 
is generated to prevent the flow of a charge current along 
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FIG. 2: (Color online) (a) DSHE and (b) ISHE signals of 
Cugg.sBio.s measured at T = 10 K. For comparison, the ISHE 
resistances of a 20 nm thick pure Cu (p = 6.3 /ifi-cm), Pt 
(10 ^if2-cm), and Cu97lr3 (14 pS7-cm) are also added in (b). 
(c) A typical AMR signal of Pyl showing the saturation of 
the magnetization above 2000 Oe for H± along the y direction 
[see Fig. 1(a)]. 



the x direction in CuBi. By inverting the probe configu- 
ration (i.e. I + <:=> V + . I~ V~ in Fig. 1(c)), one can 
also measure the direct SHE (DSHE) @, |, Il8j ; with an 
electric current in the CuBi wire, the spin accumulation 
induced at the interface between Cu and CuBi can be 
detected from the nonlocal voltage between Pyl and Cu. 

The ISHE resistance i?isHE (= Vishe/^c) a s well as 
the DSHE resistance -Rdshe for Cugg^Big 5 are plotted 
in Fig 2, with also reference signals for pure Cu, Pt, and 
Cu97ii'3. The SHE of Cu is negligibly small but, once 
only a small amount of Bi is added in Cu, the alloy shows 
a quite large SHE signal. i?isHE linearly increases with 
increasing the magnetic field and it is saturated above 
2000 Oe which is the saturation field of the magneti- 
zation, as can be seen in the anisotropic magnetoresis- 
tance (AMR) curve of Fig. 2(c). The amplitude of the 
SHE resistance is exactly the same for both the DSHE 
and ISHE, in agreement with the Onsager reciprocal re- 
lation 



18(. We note the following two points: (I)the 



sign of i?isHE for CuBi is opposite to that for Pt fj, (l, Q 
and Cu 9 7ir 3 [l8j|, and (II) the amplitude (Ai? ISH E) of 
the ISHE resistance for Cugg.sBio.s is several times larger 
than that of Pt and Cug7ir3 although their residual resis- 
tivities are almost the same (p ~ 10 /if2-cm), indicating 
that the SH angle in CuBi is much larger than in Pt and 
Culr. 

For a quantitative analysis, one has to know how much 
of the pure spin current Is generated from Pyl is ab- 
sorbed into the CuBi middle wire. As mentioned above, 
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FIG. 3: (Color online) (a) NLSV signals measured at T = 
10 K with a Cugg.sBio.s middle wire (red) and without CuBi 
wire (black). Arrows represent the magnetization directions 
of Pyl (left) and Py2 (right). Note that for the NLSV mea- 
surement the magnetic field is aligned along the easy axis (x 
direction) of the Py wires, (b) Spin diffusion length Am of 
CuBi alloys at 10 K as a function of pi mp . Closed and open 
symbols correspond respectively to the 3D and ID analyses. 



most of the spin current flowing in Cu is injected into the 
CuBi wire because of its strong SO interaction, but the 
rest of the spin current also flows toward the second Py 
wire (Py2). By measuring nonlocal spin valve (NLSV) 
signals with and without the middle wire, one can obtain 
the spin absorption rate into the CuBi wire [!, HI, EH • Fig- 
ure 3(a) shows the NLSV signals with and without the 
Cugg.sBio.s wire. With the CuBi wire, the spin accumu- 
lation signal Ai?^ ith , defined by the spin accumulation 
voltage (AVg with ) divided by I c , is reduced to ~ 50% 
compared to that without the CuBi wire (AR^ ithout ). 
The difference between ARg ithout 



and Ai?g ith depends 
strongly on the spin diffusion length of CuBi, Am, and 
can be used to extract Am, what we have done in two 
types of model: (A) the standard 1-dimensional (ID) 
model of diffusive spin transport developed by Takahashi 
and Maekawa 12 fj and already used in publications by 
some of us [si Il8l|. and (B) a 3D finite element treat- 
ment of diffusive spin transport based on an extension of 
the Valet-Fert formalism to non-collinear spins, SHE and 
3D distributions of spin and charge currents (SpinFlow 
3D [3]). More details on the two models are presented in 
supplemental material (lij ] with a list of the parameters 
which we use in the ID and 3D calculations. From the 
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data on Ai?™ hout and Ai?™ 11 we have derived A M in the 
ID and 3D models and we present the results in Fig. 3(b). 
In both models Am rapidly decays with increasing the Bi 
concentration, which expresses the contribution of the Bi 
impurities to SO-induced spin relaxation. The 3D model 
gives a longer spin diffusion length for Am > tu (thick- 



J 



ness of the CuBi wire). 

Once Am is known, we have used its respective values 
from the ID and 3D calculations to derive the SH resis- 
tivity pshe hi both approaches. In the ID approach, we 
used the following standard equation obtained by com- 
bining Eqs. (2) and (3) of Ref. [ij]: 



PSHE 



w M t M 1 - exp (-2i M /A M ) Rn {cosh(L/A N ) - 1} + 2R F {exp (L/A N ) - 1} + 2R M sinh(L/A N ) 
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where i, wm, An, and p F are respectively the distance 
between Pyl and Py2 (fixed to 1 pm), the width of the 
CuBi wire, the spin diffusion length of Cu, and the spin 
polarization of Py, while the spin resistances Rx (X = N, 
F and M) are defined in Ref. The determination of 
the shunting coefficient x has been the object of a recent 
debate [2l[ and this is one of the reasons for which we 
have introduced a 3D modeling to take automatically into 
account the shunting. As we will see later on, however, 
the larger deviation of the results of the ID model from 
those of the 3D one comes less from the shunting effects 
than from the spreading of the spin accumulation over 
the sides of the contacts in the SHE material. 

In Fig. 4(a) the SH resistivities of the CuBi alloys de- 
rived by our ID and 3D calculations are plotted as a 
function of the resistivity induced by the Bi impurities, 
Pi mp = pcuBi — Pcu- For both calculations the linear 
variation of the SH resistivity characteristic of skew scat- 
tering by dilute impurities is observed only at the lowest 
concentrations. The large deviation from linearity from 
c = 1% and above is consistent with a similar devia- 
tion from linearity in the variation of the CuBi resistiv- 
ity with the Bi concentration [see the inset of Fig. 4(a)] 
and simply reflects the departure from the dilute impu- 
rity regime. We present in Ref. fl9j scanning tunneling 
electron microscopy (STEM) and energy dispersive X-ray 
(EDX) data showing the inhomogeneous distribution on 
Bi in the concentrated alloys [22j. Note that the dilute 
impurity regime is more extended in the Culr case [lij . 
where both pshe and pimp follow linear variations up to 
Ir concentration 12% [see Fig. 4(b)]. 

From now we consider only the dilute impurity regime 
characterized by linear variations of the resistivity and 
SH resistivity with the concentration. Considering the 
slope psHE/pimp hi this regime, that is the SH angle 
an characteristic of the skew scattering by Bi impu- 
rities [18| . we remark that the 3D calculation gives a 
larger SH angle, q| d = —(0.24 ± 0.09) compared to 
a}P = -(0.12 ± 0.04). The main reason for this un- 
derestimation of the SH angle in the ID approach is 
obvious from Fig. 4(c) where one sees the 3D spread- 



ing of the spin accumulation on both side edges at the 
CuBi/Cu junction. This spreading is important when 
Am is longer than £m- Calculations with the ID model, 
in which a pure spin current is injected only vertically 
into the CuBi wire, cannot take into account the spin 
escape by lateral spreading and thus underestimate the 
SH angle. The 3D modeling allows to calculate precisely 
all the contributions coming from local distributions of 
current and shunting effects around the contacts. The 
difference between the two models is much smaller for 
the Culr alloys where Am is in general shorter than £m 
and the spin accumulation spreading is less important 
[see Figs. 4(b) and 4(d)]. The situation is similar for our 



Pt wire (Apt ~ 10 nm < tpt = 20 nm, 



, ID 



0.021 



0.024). These results are in contrast to the 



and ajf 

propositions made by Liu et al. 21( ; taking into account 
the 3D effects (shunting, spin accumulation spreading, 
etc) does not lead to a shorter Am- On the contrary Am 
is slightly longer in our 3D treatment, and thus the ID 
model underestimates the SH angle (in the present case 
by a factor of 2), not when Am is shorter than tM but 
when it is longer. Thus, the important result is the giant 
SH angle, —0.24, induced by skew scattering on Bi impu- 
rities in Cu. This SH angle is the largest value as far as 
we know, definitely above the range 0.12 — 0.15 recently 
found in Au [23| and Ta [l7|. Note that a H = -0.24, ob- 
tained by dividing pshe by Pimp, is the characteristic SH 
angle of the skew scattering on Bi. When pshe is divided 
by pcuBi(= Pcu + Pimp), the global SH angle of the alloy 
becomes —0.11. The value «h = —0.24 would be ob- 
tained by using the copper of smaller residual resistivity 
(less defect and thicker film). 

The large SH angle predicted by a recent a6-initio cal- 
culation for CuBi alloys [l4[ is consistent with our ex- 
periment but the signs are opposite. To clear up the is- 
sue, we have tried another theoretical approach adapting 
the phase shift model worked out by Fert and Levy j 1 51 ] 
for 5d impurities to the case of Bi. We start with an 
a&-initio calculation (using Quantum-ESPRESSO Q) 
of the numbers of electrons at j = 1/2 and j = 3/2 p- 
states in the conduction band on the Bi site (see Ref. [19J 
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FIG. 4: (Color online) (a) SH resistivity of CuBi alloys as a function of pi mp 
to the 3D and ID analyses 



Closed and open symbols correspond respectively 



The SH angles a| D and a™ correspond to the slopes of solid and broken lines, respectively. 



The inset shows the resistivity induced by Bi impurities pimp as a function of Bi concentration, (b) For comparison, the SH 
resistivity of Culr alloys is plotted as in (a), but the experimental data are derived from Fig. 3 in Ref. [ll|. The linear variation 
of the SH resistivity with pi mp is not limited by the solubility of Ir up to concentrations as large as 12%. (c), (d) 3D mappings 
of the spin accumulation voltage for the (c) Cugg.sBio.s and (d) Cug7iX3 SH devices calculated with SpinFlow 3D [r|. 



for more details) . By comparing with the numbers of p- 
conduction electrons on a Cu site, we find the numbers 
of electrons attracted on a Bi site and derive the j = 1/2 
and j = 3/2 scattering phase shifts (771/2 and 773/2) from 
Friedel's sum rule. The difference in these phase shifts 



comes from the SO interaction and induces the SHE. The 
phase shift 770 of the I ~ channel is also found from the 
ab-initio calculations. A strai ght forward extension of the 
formalism developed in Ref. [15J for d to p states leads 
to the following expression of the SH angle: 



-2 sin 770 {sin 771/2 sin (771 / 2 



sin 773/2 sin (773/2 



no)} 



3 (sin 7/0 + sin" 77 1/2 + 2 sin 773/2) 
I 



(2) 



The final result for the SH angle of CuBi alloy is —0.046. 
The negative sign is in agreement with our experimental 
results and the SH angle is large but not as large as in 
the experiments. However, it can be further increased, 
for example, by taking into account electron correlation 
effects which may significantly enhance the SO splitting 
as recently proposed by Gu et al [23I ] . Indeed, we have 
checked that an increase of the splitting between the 771/2 
and 773/2 phase shifts by a factor of about 3 keeps the right 
sign but would bring the SH angle beyond —0.26. 

Finally one wonders if interfacial effects could also con- 
tribute to our experimental results. The possible compe- 
tition between the SHE and the interfacial Rashba effects 



has indeed been debated recently [25j, |26j for the inter- 
pretation of experiments with 2 ~ 3 nm thick Pt layers. 
For our thicker CuBi films (20 nm), the SHE contribu- 
tion can be expected to be largely predominant. Also 
some enhancement of the SHE at interfaces 27| can be 
neglected in a 20 nm thick metallic film in which the 
interface states cannot exceed a couple of atomic layers. 

In conclusion, we find that a small amount of Bi im- 
purities (< 0.5%) in Cu induces a large SHE. We have 
analyzed our data using both the classical ID model of 
previous SHE studies and a 3D finite element treatment 
of spin transport. A large SH angle is derived from both 
models. It is definitely larger in the 3D model (—0.24). 
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Such a difference between ID and 3D models is not sur- 
prising since the 3D model can treat more accurately the 
unavoidable approximations of the ID model. Harness- 
ing such a giant SHE to produce or detect spin currents 
will be probably more and more used in novel generations 
of spintronic devices not necessarily based on magnetic 
materials. 
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